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We study gravitational waves in the black string Randall-Sundrum braneworld. We present 
a reasonably self-contained and complete derivation of the equations governing the evolution of 
gravitational perturbations in the presence of a brane localized source, and then specialize to the 
case of spherical radiation from a pointlike body in orbit around the black string. We solve for the 
resulting gravitational waveform numerically for a number of different orbital parameters. 



I. INTRODUCTION 

The Randall-Sundrum (RS) braneworld model postu- 
lates that our universe is a 4-dimensional hypersurface 
embedded in 5-dimensional space [H I2 E]- One of the 
most remarkable features of the model is that in the 10 
years since its introduction, no one has found any ev- 
idence of tension with current observations or tests of 
gravitational phenomena. This is despite the fact that 
the model incorporates a large extra dimension through 
which gravity is allowed to propagate. The RS model's 
viability stems from the fact that it alters conventional 
general relativity (GR) on scales smaller than the curva- 
ture scale of the bulk spacetime i. Hence, if one chooses (. 
to be sufficiently small the RS model is indistinguishable 
from GR in many experimental or observational situa- 
tions. 

The principal virtue of the RS scenario is also a bit of a 
detriment: In order to constrain or refute the model one 
has to look at increasingly smaller scale phenomena. The 
most direct test comes from laboratory measurements of 
the gravitational force between two masses, which yields 
^ < 50 /xm lU [S] . One can also examine high energy cos- 
mological phenomena to derive observable consequences 
of the braneworld paradigm. The idea is that when the 
Hubble radius becomes smaller than the bulk curvature 
Ht < 1, the physical scale of all interesting gravitational 
interactions are also smaller than £. So in these epochs, 
one expects the RS corrections to GR to become domi- 
nant. The spectrum of tensor perturbations in the high 
energy radiation RS era has been calculated and shown 
to the be consistent with the GR result (with minor mod- 
ifications) [6l[7l|8]. On the other hand, the spectrum of 
scalar density perturbations is found to be enhanced over 
the GR expectation in the early universe, which could 
lead to an overproduction of primordial black holes in 
the RS model The behaviour of scalar perturba- 
tions during inflation has also been considered, and it 
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the same shadow 




FIG. 1: A schematic of the black string braneworld. The 
black string is a line singularity extending from our brane to 
a shadow brane in the bulk. There is a periodic identification 
such that the background geometry is symmetric about our 
brane. The singularity is covered by a warped horizon, and 
the brane separation must be large enough to be consistent 
with solar system tests of GR yet small enough to avoid the 
Gregory-Laflamme instability [121 1131 [T^ . In this paper, we 
consider the situation where the string is perturbed by a small 
orbiting body on one of the branes, here depicted as a small 
black hole. 



was found that there were very small corrections to the 
power spectrum of primordial fluctuations |10l lllj . 

Another possible means of constraining the RS model 
is with gravitational waves (GWs) with wavelengths A < 
£. These naturally probe gravitational interactions in 
the regime where RS effects should be important, and 
can be viewed as the dynamical counterpart to the static 
laboratory tests of Newton's Law mentioned above. It 
is useful to have a concrete model of the generation and 
propagation of these short-wavelength GWs in order to 
determine if they have sufficient amplitude to be observed 
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by real GW detectors. To that end, we have previously 
considered the behaviour of gravitational perturbations 
around a black string braneworld (see Fig.[T]) and showed 
that a generic feature of the signal involves a long-lived 
oscillatory tail composes of a discrete spectrum of modes 
whose wavelengths are less than I [14 . The amplitude 
of spherical radiation emitted by a black string being 
perturbed by an orbiting small body was estimated in 
[T5] . This type of GW radiation is considered to be a 
possible source for very-high frequency GW detectors [T6| 

nzmH]. 

Our purpose in this paper is to present the details 
of the black string perturbative formalism utilized in 
the Letters [Ml [15]. This is the subject of jP-fvn 



We also describe how to numerically calculate the GWs 
sourced by a "point particle" orbiting the black string, 
and present the results of a number of simulations in 

^vng 



II. A GENERALIZED RANDALL-SUNDRUM 
TWO BRANE MODEL 

In this section, we present a generalized version of the 
Randall-Sundrum two brane model in a coordinate in- 
variant formalism. Our treatment represents a general- 
ization of the work of Shiromizu et al. ^19^. We begin by 
outlining the geometry of the model, the action governing 
the dynamics, and the ensuing field equations. We then 
specialize to the black string braneworld model, which 
will be perturbed in the next section. 



A. Geometrical framework and notation 

Consider a (4+l)-dimcnsional manifold (A^,g), which 
we refer to as the 'bulk'. One of the spatial dimensions 
of M. is assumed to be compact; i.e., the 5-dimensional 
topology is X S. We place coordinates on M. so 
that the 5-dimensional line element reads: 



(isg = gAsdx dx 



(2.1) 



We assume that there is a scalar function <f> that uniquely 
maps points in M. into the interval / = (— d, Here, 
d is a constant parameter that is one of the fundamen- 
tal length scales of the problem. The gradient of this 
mapping 9^<i> satisfies 



> 0, 



(2.2) 



and is tangent to the compact dimension of Ai. This 
scalar function defines a family of timelike hypersurfaces 
<i>(x"^) — Y, which we denote by Sy. The two submani- 
folds at the endpoints of /, and S-d, are periodically 
identified. 

Let us now place 4-dimcnsional coordinates z" on each 
of the Sy hypersurfaces. These coordinates will be re- 
lated to their 5-dimensional counterparts by parametric 



equations of the form: 
following basis vectors 

A _ dx^ 
riAe^ = 



x^{z"). We then define the 



-1. 



n riA 



(2.3) 

The tetrad is everywhere tangent to Sy, while is 
everywhere normal to Sy. The projection tensor onto 
the Ey hypersurfaces is given by 

qAB ^ QAB -riAriB, n'^qAB^O- (2.4) 

From this, it follows that the intrinsic line element on 
each of the Ey hypersurfaces is 

dsl = qapdz"^ dz^ , = e^e^qAB = e^e^gAB- (2.5) 

The object behaves as a tensor under 4-dimcnsional 
coordinate transformations z" z"{z^) and is the in- 
duced metric on the Ey hypersurfaces. It has an inverse 
q"'' that can be used to define e'^: 



^A 



(2.6) 



Generally speaking, we define the projection of any 5- 
tensor Tab onto the Sy hypersurfaces as 



Tap — e^epTAB, 



(2.7) 



where the generalization to tensors of other ranks is ob- 
vious. The 4-dimensional intrinsic covariant derivative of 
Taj3 is related to the 5-dimensional covariant derivative 
of Tab by 



[^M, = ete'^e^^Aqftq^TBC: 



(2.8) 



where the notation [■ ■ ■\q means that the quantity inside 
the square brackets is calculated with the qap metric. 

Finally, the extrinsic curvature of each Sy hypersur- 
face is: 

Kab ^ q^ycriB ^ \£,iqAB ^ Kba, ti^Kab^O, 

(2.9) 



B. The action and field equations 

We label the hypersurfaces aX Y = y+ — Q and Y = 
y_ = +d as the 'visible brane' and 'shadow brane' 
E", respectively. Our observable universe is supposed to 
reside on the visible brane. These hypersurfaces divide 
the bulk into two halves: the lefthand portion A^l which 
has y G (— d, 0), and the righthand portion which has 
y G {0,+d). The action for our model is: 



S 



1 

2^ 



Ml 



Ml 



(5)i?-2A, 



1 
2 

Mb 



1 

24 



(5)i?-2A, 



Mb 



(2.10) 
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In this expression, Kg is the 5-dimensional gravity matter 
coupling, A5 = —6k^ is the bulk cosmological constant, 
A* — ibG/c/Kg are the brane tensions, and i = 1/fc is the 
curvature length scale of the bulk. Also, is the La- 
grangian density of matter residing on S*, while >Cl and 
Cn are the Lagrangian densities of matter living in the 
bulk. Note that the visible brane in our model has posi- 
tive tension while the shadow brane has negative tension. 

The quantity [K]^ is the jump in the trace of the 
extrinsic curvature of the Sy hypersurfaces across each 
brane. To clarify, suppose that dAi^ and are the 

boundaries of and Mu coinciding with S , respec- 
tively. Then, 



[K]+ = q'-^K^p 



dM- 



(2.11a) 
(2.11b) 



We can now write down the field equations for our 
model. Setting the variation of S with respect to the 
bulk metric g^^ equal to zero yields that: 



Gab ~ Gk'gAB = 4 [e{+y)T% + e{~y)T\s] 
2 Si^Ci^^n) 



/^-^L,R. 



AB 



(2.12) 



Meanwhile, variation of S with respect to the induced 
metric on each boundary yields 

Q±s = {[Kab] ± 2kqAB + 4{Tab - \TqAB)]^ = 0, 



^AB 



5q' 



(2.13) 



Here, the {•••}* notation means that everything in- 
side t he curly brackets is evaluated at We see that 
(2.121 are the bulk field equations to be satisfied by the 
5-dimensional metric Qabi while (2.131 are the bound- 



ary conditions that must be enforced at the position of 
each brane. Of course, (2.131 are simply the Israel junc- 
tion conditions for thin shells in general relativity. In 
a braneworld context, the symmetric versions of these 
equations first appeared in Shiromizu et al. |19j . 

In what sense is our model a generalization of the RS 
setup? The original Randall-Sundrum model exhibited a 
Z2 symmetry, which implied that A^l is the mirror image 
of A^R. Also, in the RS model the bulk was explicitly 
empty. However, since we allow for an asymmetric dis- 
tribution of matter in the bulk, we explicitly violate the 
Z2 symmetry and bulk vacuum assumption. 



C. The black string braneworld 

We now introduce the black string braneworld, which 
is a Z2 symmetric solution of (2.121 and (2.131 with no 
matter sources: 



R 



0. 



(2.14) 



Here, we use = to indicate equalities that only hold in 
the black string background. The bulk geometry for this 
solution is given by: 



fir. 



-f{r)df 
= 1 



1 



fir) 
2GM/r, a{y) 



-k\v\ 



(2.15) 



Here, M is the mass parameter of the black string and 
G = lp\/M-p\ is the ordinary 4-dimensional Newton's con- 
stant. The function <!> used to locate the branes is trivial 
in this background: 



(2.16) 



which means that the branes are located at y = and 
y = d, respectively. The Ey = T.y hypersurfaces have 
the geometry of Schwarzschild black holes, and there is 
5-dimensional line- like curvature singularity at r ~ 0: 



nABCD p 

-K HabCD 



48G2M2e4'=l^l 



40fc^ 



(2.17) 



Note that the other singularities at ?/ = ±cx3 are excised 
from our model by the restriction y = Y ^ (— d, d], so we 
will not consider them further. 

Finally, note that the normal and extrinsic curvature 
associated with the Sy hypersurfaces satisfy the follow- 
ing convenient properties: 



n-A = dAy: 



h^Vati^ 



0, K 



AB 



-kq_ 



AB- 



(2.18) 



These expressions are used liberally below to simplify 
formulae evaluated in the black string background. 



III. LINEAR PERTURBATIONS 

We now turn attention to perturbations of the black 
sting braneworld. Our treatment will be a reformulated 
and generalized version of the original Randall-Sundrum 
work [1, 2, and the seminal contribution of Garriga and 
Tanaka EUl. 



A. Perturbative variables 

We are ultimately interested in the behaviour of gravi- 
tational waves in this model, which are described by fluc- 
tuations of the bulk metric: 



9AB gAB + hAB, 



(3.1) 



where Hab is understood to be a 'small' quantity. The 
projection of Hab onto the visible brane is the observable 
that can potentially be measured in gravitational wave 
detectors. But it is not sufficient to consider fluctuations 
in the bulk metric alone — to get a complete picture, we 
must also allow for the perturbation of the matter content 
of the model as well as the positions of the branes. 
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Obviously, matter perturbations are simply described 
by the T^^, T^g, and T^g stress-energy tensors, which 
are considered to be small quantities of the same order as 
hAB- On the other hand, we describe fluctuations in the 
branc positions via a perturbation of the scalar function 



(3.2) 



Here, ^ is a small spacetime scalar. Recall that the po- 
sition of each brane is implicitly defined by ^{x^) = y±. 
Hence, the brane locations after perturbation are given 
by the solution of the following for y: 



v=y± 



{y - y±)dyi 



y=y± 



y±- 



(3.3) 



However, note that y — y± is of the same order as ^, so at 
the linear level the new brane positions are simply given 

by 



y = y±-^ 



(3.4) 



Hence, the perturbed brane positions are given by the 
brane bending scalars: 



0. 



(3.5) 



Note that because and ^~ are explicitly evaluated at 
the brane positions, they are essentially 4-dimensional 
scalars that exhibit no dependence on the extra dimen- 
sion. 

Having now delineated a set of variables that parame- 
terize the fluctuations of the black string braneworld, we 
now need to determine their equations of motion. 



B. Linearizing the bulk field equations 



First, we linearize the bulk field equations (2.12 1 about 
the black string solution. Notice that (2.12) only de- 



pends on the bulk metric and the bulk matter distri- 
bution. Hence, the linearized field equations will only 
involve Hab , and . The actual derivation of the 
equation proceeds in the same manner as in 4-dimcnsions, 
and we just quote the result: 



y'^^chAB - y^^AhBC - V^Vs/lAC 



'hAi 



-2Kl^\t, (3.6) 



where 



^AB 



e{+y){T%-lT^gAB) 



ei-y)iT^B-lT^9AB). (3.7) 



The wave equation (3.6 1 is valid for arbitrary choices of 



gauge and generic matter sources. If we specialize to the 
Randall-Sundrum gauge 



^""hAi 



0, 



0, fiAB = ele^B^ap, (3.8) 



eq. (3.6 1 reduces to 



kAB^'^hcD + {GMaf{£l - 4k^)hAB 

= -2{GMa)'4E\t, (3.9) 
where we have defined the operator 



A 



- 2(4)i? 



C D^ 
A B \ 



= (GM) 



^A^B 



2Raf3 



ezej. (3.10) 



Here, ^^^Racbd is the Riemann tensor on Sy, which can 
be related to the 5-dimcnsional curvature tensor via the 
Gauss equation 

^'^^RmNPQ = qAigNQplQ^ABCD + 2Km[pKq^n. (3.11) 



On the second line of ( 3.10 ) the 4-tensor inside the square 
brackets is calculated using q^fj. We can re-express this 
object in terms of the ordinary Schwarzschild metric gap, 
which is conformally related to qap via the warp factor: 



9al3 



dz"dz>^ 



qap 

-fdt^ + f 



(3.12a) 
(3.12b) 



The quantity in square brackets on the third line of ( 3. 10 1 



is calculated from ga/s-^ One can easily confirm that 
Aab'^^ is 'y-independent' in the sense that it commutes 



[(4)A, 



AB , £n] — 0. 



(3.13) 



In addition, the (GM)^ prefactor makes Aab'~''^ dimen- 
sionless. 



Notice that the lefthand side of (3.9) is both trace- 
less and manifestly orthogonal to n^, which implies the 
following constraints on the bulk matter: 



^AB — ^A^B^aP ; 1 

In other words, our gauge choice is inconsistent with 
bulk matter that violates these conditions. If we wish 
to consider more general bulk matter, we cannot use the 
Randall-Sundrum gauge. 



Q/3 y^bulk 



= 0. 



(3.14) 



C. Linearizing the junction conditions 

Next, we consider the perturbation of the junction con- 
ditions (2.13). These can be re- written as 



Qab = [ik^iAnB) - "-(A|»^'^Vcn.|B)] 

±kqAB + 4{TAB-lTqAB)}^ =0. (3.15) 



^ Unless otherwise indicated, for the rest of the paper any tenso- 
rial expression with Greek indices should be evaluated using the 
Schwarzschild metric Qcfj- 
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We require that Q^g vanish before and after perturba- 
tion, so we need to enforce that the first order variation 
^Qab equal to zero. 

In order to calculate this variation, we can regard the 
tensors Q^g as functionals the brane positions (as de- 
fined by <&), the brane normals ua, the bulk metric, and 
the brane matter: 

QaB = QABi'^^'^M,gMN,T^jN), (3.16) 

from which it follows that 



„ 5^ + one 

[ (5$ 6nc 

, SQab . , SQab .rj, \^ 

+ ^ (>gcD + ^ — 0-lcD> ■ 

dgcD 01 CD Jo 



(3.17) 



The {• • • }J notation is meant to remind us that after 
we have calculated the variational derivatives, we must 
evaluate the expression in the background geometry at 
the unperturbed positions of the brane. 

We now consider each term in (3.171. For simplic- 



ity, we temporarily focus on the positive tension visible 
brane and drop the -I- superscript. The first term rep- 
resents the variation of Qab with brane position, which 
is covariantly given by the Lie derivative in the normal 
direction: 



Ms 



(5$ 



6<i>\ ={-C£nQAB}o 



(3.18) 



But the Lie derivative of Qab vanishes identically in the 
background geometry, so this term is equal to zero. 
The second term in (3.171 represents the variation of 



Qab with respect to the normal vector. Making note of 
the definition (2.3 1 of in terms of <&, as well as (5<i> = ^ 
I'^VyiC 0; arrive at 



and 



Sua = Va'?, n Sua — 0. 



(3.19) 



Notice that since the normal itself must be continuous 
across the brane, we have [6nA] = 0. After some algebra, 
we find that the variation of the junction conditions with 
respect to the brane normal is non-zero and given by 



^<5nc}^ = 2g2gfVcV^S. (3.20) 



The third term in (3.171 is the variation with the bulk 



metric itself 6gAB ~ It-ab- Calculating this is straight- 
forward, and the result is: 



^-^Sgco] =\[£nhAB]+'ikhAB- (3.21) 

ogcD J 



The last variation we must consider is with respect to the 
brane matter fields, which is trivial: 

^^^^""-STcd] =nl{TAB-lTqAB). (3.22) 



5T, 



CD 



So, we have the final result that 

5Q% = {2q'iq^VcVD£. + \[£nhAB] 

±2khAB + 4 {Tab - ITqab) }o - 0. (3.23) 



If we take the trace of SQ^g — 0, we obtain 

„AB-r-, Y7 f± 1 ^2rp± 



(3.24) 



These are the equations of motion for the brane bending 
degrees of freedom in our model, which are seen to be 
directly sourced by the matter fields on each brane. 



D. Converting the boundary conditions into 
distributional sources 



We can incorporate the boundary conditions 5Q^ 



AB 



directly into the Hab equation of motion as delta- 
function sources. This is possible because the jump in 
the normal derivative of Hab appears explicitly in the 
perturbed junction conditions. This procedure gives 



^AB^^hcD — fl^hAB 



= -2{GMafKl 
Here, we have defined 



vbulk 



AB 



AB 



(3.25) 



i2 = -(GMa)^ 

Ei f T'^ 1 T^i 



2k 



^^A^5(y-2/,)-4fc2 



QAB 



TQlqE^c^D^^. (3.26) 



If we integrate the wave equation (3.251 over a small re 



gion traversing either brane, we recover the boundary 
conditions (3.23 1. 



Together with the gauge conditions 

iB — q' ' AIT-CB 



n^hAB = q^^'^AhcB = = q^^hAB, 



(3.27) 



(3.241 and (3.25) are the equations governing the pertur- 



bations of our model. 

IV. KALUZA-KLEIN MODE FUNCTIONS 

A. Separation of variables 

As mentioned above, we have that 

[^ab'~^^ , £n]hcD — 0; 



(4.1) 



i.e., Aab'^^ is independent of y when evaluated in the 
{t,r,6,(j),y) coordinates. This suggests that we seek a 
solution for Hab of the form 



Hab = ZHab, fi^Z — [I Z, 



(4.2) 
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where. 



Case 1: light modes 



= £nhAB and = q^VAZ; 



(4.3) 



that is, Z is an eigenfunction of ji^ with eigenvalue /i^. 
The existence of the delta functions in the jj? operator 
means that we need to treat the even and odd parity 
solutions of this eigenvalue problem separately. 



B. Even parity eigenfunctions 

If Z{—y) — Z{y), we see that Z satisfies the following 
equations in the interval y € [0, d]: 

m'Z{y) ^ -a^{y){dl - ie)Z{y), 

0^[{dy + 2k)Z{y)]±, (4.4) 
fi = GMm. 

There is a discrete spectrum of solutions to this eigen- 
value problem that are labeled by the positive integers 
71= 1,2,3...: 

Zniy) = a~'[Yi{mn£)J2{mn£e''\y\) 

- MmJ)Y2{m^£e''^y% (4.5) 

where is a constant, and m„ = Hn/GAI is the n**^ 
solution of 

yi(m„^)Ji(TO„fe*'-'') = Ji{mne)Yi{m„ee''^). (4.6) 

There is also a solution corresponding to toq — — 0, 
which is known as the zero- mode: 



Zo{y) = ^e" 



■2k\v\ 



"0 



-2kdy/2 



(4.7) 



Hence, there exists a discrete set of solutions for bulk 
metric perturbations of the form h^^^ = Z„{y)h^^g{z°'). 
When n > these are called the Kaluza-Klein (KK) 
modes of the modes, and the mass of any given mode is 
given by the m„ eigenvalue. The a„ constants are deter- 
mined from demanding that {Z„} forms an orthonormal 
set 

Srnn = / dy a '^{y)Z,n{y)Zn{y). (4.8) 
These basis functions then satisfy: 

oo 

S{y ~y±)^J2 a-'Z^{y)Z^{y±). (4.9) 

n=0 

This identity is crucial to the model — inspection of 
(3.251 reveals that the brane stress energy tensors ap- 



pearing on the righthand side are multiplied by one of 
S{y — y±). Hence, brane matter only couples to the even 
parity eigenmodes of ft^. 



It is useful to have simple approximate forms of the 
Kaluza-Klein masses and normalization constants. These 
are straightforward to derive for modes that are 'light' 
compared to mass scale set by the AdS^ length parame- 
ter: 

m„£<l. (4.10) 
Let us define a set of dimensionless numbers a;„ by: 

(4.11) 



kd 



Then for the light modes, we find that Xn is the n*'' zero 
of the first-order Bessel function: 



Ji(x„) = 0. 



(4.12) 



Also for light modes, the normalization constants reduce 
to 

a„^2Vie^''^\Joixn)\/nx„, n>0. (4.13) 

Actually, it is more helpful to know the value of the KK 
mode functions at the position of each brane. We can 
parameterize these as 

^n(2;±) = Vfce-'='^z±, n>0. (4.14) 

For the light Kaluza-Klein modes, the dimensionless 
are given by 



\Joi^n)\'' 



(4.15) 



Case 2: heavy modes 

At the other end of the spectrum, we have the heavy 
Kaluza-Klein modes 



runi » 1. 

Under this assumption, we find 
riTT 



'kd ' 



Zniy) 



c 

1 

^ — kd 



gfcls/l _ 1 



mr- 



,kd 



1 



(4.16) 

(4.17a) 
(4.17b) 
(4.17c) 



(Strictly speaking, an asymptotic analysis leads to formu- 
lae with n replaced by another integer n' on the righthand 
sides of Eqns. (4.171. However, we note that for even 



parity modes, n counts the number of zeroes of Zn{y) in 
the interval y S (0, d), which allows us to deduce that 
n' = 71.) Unlike the analogous quantities for the light 
modes, shows an explicit dependence on the dimen- 
sionless brane separation d/£. 
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C. Odd parity eigenfunctions 

As mentioned above, brane matter only couples to 
Kaluza-Klein modes with even parity. But a complete 
perturbative description must include the odd parity 
modes as well; for example, if we have matter in the 
bulk distributed asymmetrically with respect to y = 
(i.e. T^g ^ T^b) modes of either parity will be excited. 
Hence, for the sake of completeness, we list a few prop- 
erties of the odd parity Kaluza-Klein modes here. 

Assuming Z(—y) = —Z(y), we have: 



m'Z{y)^-a^y)idl-Ak')Z{y), 
= Z(y+) = Z(y_). 



(4.18) 



Again, we have a discrete spectrum of solutions, this time 
labeled by half integers: 

Zn+i{y) = a;;i[r2(m„+i£)J2(m„+ife'=l^l) 

- J2(m„+i^)y2(m„+i^e^-|^l)]. (4.19) 
The mass eigenvalues are now the solutions of 

r2(m„+i^)J2(m„+i&^-'^) = 

Mm^+ii)Y,{m^+.£e''''). (4.20) 

Proceeding as before, we define 



For light modes with to„^i^ <^ 1, x,- 
of the second-order Bessel function: 



(4.21) 



is the n*^ zero 



(4.22) 



Taken together, (4.12) and (4.221 imply the following for 
the light modes: 



"^l < m3/2 < ma < TOg/j < • • • 



(4.23) 



i.e., the first odd mode is heavier than the first even 
mode, etc. 

Finally, we note that since the odd modes vanish at 
the background position of the visible brane, it is impos- 
sible for us to observe them directly within the context 
of linear theory. This can change at second order, since 
brane bending can allow us to directly sample regions of 
the bulk where Z^^i ^ 0. However, this phenomenon is 
clearly beyond the scope of this paper. 



D. Stability criterion 

Finally, as discussed in detail elsewhere [14], the black 
string braneworld will be perturbatively stable if the 
smallest KK mass satisfies 



Under the approximation that the first mode is light 
{x\e~^'^ ^ 1) and using G = €pi/Mpi, this gives a re- 
striction on the black string mass 



or equivalently, 
M 



— > 8 X 10"^ , „ , 
Mq ^ V 0.1 mm 



(4.25) 



(4.26) 



If we take ^ = 0.1 mm, then we see that all solar mass 
black holes will in actuality be stable black strings pro- 
vided that d/i < 19. 



V. RECOVERING 4-DIMENSIONAL GRAVITY 

Let us now describe the limit in which we recover gen- 
eral relativity. (Garriga and Tanaka [5D] first considered 
this problem in Minkowski space, but the approach em- 
ployed here is somewhat different.) We assume there are 
no matter perturbations in the bulk and on the hidden 
brane; hence, we may consistently neglect the odd par- 
ity Kaluza-Klein m odes. By virtue of the brane bending 
equation of motion (3.24), we can consistently set ^~ = 0. 



tion in front of E 



Furthermore, (4.9) can be used to repla ce the delta func 



AB 



in equation (3.25). We obtain 



^AB^^hcD ~ P?hAB 

oo 

= -2{GMfnl->:\^ Zn{y+)Zn{y). (5.1) 

n=0 

We now note that for e^^'^ <^ 1, 
Zo(t/+)-\/fc(l-e-2'='^)-i/2»Z„(y+), n>0. (5.2) 

That is, the n > terms in the sum are much smaller 
than the 0*^ order contribution. This motivates an ap- 
proximation where the n > terms on the righthand side 
of (5.1) are neglected, which is the so-called 'zero- mode 



truncation'. 

When this approximation is enforced, we find that hj^B 
must be proportional to Zo^y); i.e., there is no contribu- 
tion to hj^B from any of the KK modes. Hence, we have 
fi^hAB — 0. The resulting expression has trivial y depen- 
dence, so we can freely set y = j/+ to obtain the equation 
of motion for Hab at the unperturbed position of the vis- 
ible brane: 



^AB n-CD 



-2{GM)\li:+sZ^,{y, 



(5.3) 



But we are not really interested in h'^^, the physically 
relevant quantity is the perturbation of the induced met- 
ric on the perturbed brane, which is defined as the vari- 
ation of 



III = GMmi > fic^ 0.4301. 



(4.24) 



Iab = [9AB - nAUBY 



(5.4) 
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We calculate Sq\g in the same way as we calculated 
^Qab above (except for the fact that qab shows no ex- 
phcit dependence on T^^)- 

= 1^^*+ ^'^^^"/o ■ 

These variations are straightforward, and we obtain: 

^IaB = ^AB = ^AB + 2fc^+glB 

~ [UA^ B + ne"^ A)t , (5.6) 

where all quantities on the right are evaluated in the 
background and at the unperturbed position of the brane. 
Note that JiABn^ 7^ 0, which reflects the fact that ua is 
no longer the normal to the brane after perturbation. 
We now define the 4-tensors 



h+ - p'^p^h+ 



Here, h'^a is the actual metric perturbation on the visible 
brane. Note that this perturbation is neither transverse 
or tracefree: 



(5., 



We can now re-express the equati on o f motion (5.3) in 
terms of h'^^ instead of /ij^^ using (5.6 1. Dropping the + 
superscripts, we obtain 



2ZIkI 



1 



aP 



1 



Pi 
k 



fdaP 



3 V ■ 2^2 

{6k-AZl)V^Vp^, (5.9) 



where we have defined 



(5.10) 



In obtaining this expression, we have made use of the ^ 
equation of motion: 



(5.11) 



Note that we still have the freedom to make a gauge 
transformation on the brane that involves an arbitrary 4- 
dimensional coordinate transformation generated by 7]^: 



l-ap 



T-aP 



^aVp + ^pVa- 



(5.12) 



We can use this gauge freedom to impose the condition 
Vph'^a - ^y^h^p = i2Zl - 3fc)VaC. (5.13) 

Then, the equation of motion for 4- metric fiuctuations 
reads 



a. p '''7(5 



167rG 



Tap — 



1 + ^BD 

3 + 2wbd 



T\ga^ 



(5.14) 



where we have identified 



-Bo=2(e^'^/'-l)> 



(5.15) 



We see that (5.141 matches the equation governing grav- 



itational waves in a Brans-Dicke theory with parameter 
t^BD. Hence in the zero- mode truncation, the perturba- 
tions of the black string braneworld are indistinguishable 
from a 4-dimensional scalar tensor theory. 

Note that (5.141 must hold everywhere in our model. 



so we can consider the situation where our solar system is 
the perturbative brane matter located somewhere in the 
extreme far-field region of the black string. The forces 
between the various celestial bodies will be governed by 
(5.141 in the Rap-yS ~ limit. In this scenario, solar 
system tests of general relativity [2T place bounds on 
the Brans-Dicke parameter, and hence d/t. 



> 4 X 10"* 



d/l > 5. 



(5.16) 



This lower bound on the dimensionless brane separation 
will be an important factor in the discussion below. 



VI. SPHERICAL WAVES ON THE BRANE 

In this section, we specialize to the situation where 
there is perturbative matter located on one of the branes 
and no other sources. Unlike Sec. |V] our interest here is 
to predict deviations from general relativity, so we will 
not use the zero- mode truncation. Principally for rea- 
sons of simplicity, we will focus on spherically symmetric 
radiation, which is a channel unavailable in the standard 
4-dimensional setup. 



A. Mode decomposition 



To begin, we make the assumptions 



E'Xt = 0, and = or = 0; 



(6.1) 



i.e., we set the matter perturbation in the bulk and one of 
the branes equal to zero. Note that due to the linearity 
of the problem we can always add up solutions corre- 
sponding to different types of sources; hence, if we had 
a physical situation with many different types of matter, 
it would be acceptable to solve for the radiation pattern 
induced by each source separately and then sum the re- 
sults. 

We decompose Hab as 



hAi 



^ e^e^ 2^Z„{y)Zn{y±)h)^f^. 



(6.2) 



Here, C is a normalization constant (to be specified later) 
with dimensions of (mass)"**, and the expansion coeffi- 
cients /i^'^^ are dimensionless. We define a dimensionless 
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brane stress-energy tensors and brane bending scalars by 



(6.3) 



Omitting the ± superscripts, we find that the equation 
of motion for h^'^jj is 



{GMf 



- 2 (e„^ - lega,„) - 4(GM)2V„V0|, (6.4) 
while the equation of motion for ^ is 

V"V„|=ie. (6.5) 
We also have the conditions 



'q/3 ■ 



(6.6) 



Note that in all of these equations, all 4-dimensional 
quantities are to be calculated with the Schwarzschild 
metric g^^. In particular, Q = g^^Oap- 



B. The radiative s-wave channel 

We now turn our attention to solving the coupled sys- 
tem of equations (6.4| and (6.51 for ageneric source Qafi- 
The symmetry of the background geometry dictates that 
we decompose the problem in terms of spherical harmon- 
ics: 



^ /A 

, (n) _ "a/3 
"'^ " V4^ 



(6.7a) 



l—l rn— — l 
oo I 10 

E EK™]"/3^^'"'(6-7b) 



0Q/3 — 



e 



a/3 

47r 



oo ; 10 



E E Et^.™]"^©^^- 



/^l m——l i—1 



Here, [^;|^^]c(/3 are the tensorial spherical harmonics in 4 
dimensions. The terms with ^ > in this decomposition 
can be quite involved, so for the purposes of this paper 
we concentration on the spherically symmetric s-wave 
(Z = 0) sector, which is described by ^^ai3^\ ^'^d 



contribution to the metric pertur- 



We write the 
bation as 

^(n^s) ^ |_|^ ^^^^ _ 2|_|^ t(^r^) + Hg r^ra -l- K 7^/3, (6.8) 

where we have defined the orthonormal vectors 

t"-r'/'a„ r" = /i/2a„ (6.9) 

which are pointing in the time and radial directions, re- 
spectively; and 

lap = gal3 + tatf3 - rarfj, t^JafS = r°'7aP ^ 0. (6.10) 



which is the induced metric on the 2-spheres of constant 
r and t. Each of the expansion coefficients is a function 
of t and r; i.e., — Hi{t, r) and K — K{t, r). Notice that 

the condition that '^h'"^} is tracefree implies 



K= i(Hi-H3). 



(6.11) 



Before going further, it is useful to define dimensionless 
coordinates: 



P 



GM' 



GM' 



2; = p-f 21n(^-l) . (6.12) 



Then, when our decompositions ( |6.7[ ) are substituted into 
the equations of motion, we find that all components of 
the metric perturbation are governed by master variables 



2p3 



2 + fi^p^ 



ip = P- 



(6.13) 



Both ip — iP{t, x) and tp ~ (p(r, x) satisfy simple wave 
equations: 

(9^ -dl + V,p)^P = S^+i^, (6.14a) 
{d^,-dl + V^)^ = S^. (6.14b) 

The potential and matter source term in the equation 
are: 



f 



p3(2 + p3^2)2 

2fp' 
3(2 + /xV)' 



,4„6 



~2Ap^p^ + 36p^p^ 
p(2 + ^V)5r(2Ai-H3A3) 



-f 6(/iV-4)/A2 



(6.15a) 



(6.15b) 



Here, we have defined the following three scalars derived 

(s) 

from the dimensionless stress-energy tensor S^,^: 

A2 = -eiit°'r'^ 



A3 - 



(6.16a) 
(6.16b) 
(6.16c) 



The potential and source terms in the brane-bending 
equation are somewhat less involved: 



Pf 



drAi. 



(6.17) 



Finally, the interaction operator is 



I = 



(2 + /x2p3)' 



[6/p25p + (^V-6p+8)] . (6.18) 
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Inversion formulae 



Assuming that we can solve the wave equations (6.14) 



for a given source, we need formulae that allow us to 
express H^, K in terms of ip and (p in order to make 
gravitational wave prediction. This can be derived by 



inverting the master variable definitions (6.13) with the 



aid (6.14). The general formulae are actually very com- 



plicated and not particularly enlightening, so we do not 
reproduce them here. Ultimately, to make observational 
predictions it is sufficient to know the form of the metric 
perturbation far away from the black string and the mat- 
ter sources, so we evaluate the general inversion formulae 
in the limit of p — > oo and with = 0: 



H2 



do 



do 



1 



2 1 i'P 
MP V P 



(6.19) 



Note that these do not actually complete the inversion; 
in most cases, a quadrature is also required to arrive at 
the final form of the metric perturbation. 



particle follows a geodesic along the brane: 



dz" dz^ 



^'^ ^ dr] di] ' '^"^ drj drj ^ ^ ^ 



where r^^[g^] are the Christoffel symbols defined with 
respect to the metric. 



We note that the above formulae make explicit use of 
the induced brane metrics q^^. However, all of our per- 
turbative formalism is in terms of the Schwarzschild met- 



ric gai3, especially the definition of the A^ scalars (6.16) 



Hence, it is useful to translate the above expressions us- 
ing the following definitions: 



dz^ 
~dX' 



-1 = gapU^W 



(7.5) 



Then, the stress-energy tensor and particle equation of 
motion become 



Mp 
a± 



(54(z^ 



=^^UaUpdX, u"Vq,W^ = 0. 



'9 



(7.6) 

Note that the only difference between the stress-energy 
tensors on the positive and negative tension branes is an 
overall division by the warp factor. 

By switching over to dimensionless coordinates, trans- 
forming the integration variable to r from A, and making 
use of the spherical harmonic completeness relationship, 
we obtain 



VII. POINT PARTICLE SOURCES ON THE 
BRANE 

We now specialize to the situation where the perturb- 
ing brane matter is a "point particle" located on one of 
the branes. We take the particle Lagrangian density to 
be 




S'^{z^' - z^) dz^ dz^ 
\/-~q drj drj 



(7.1) 



In this expression, 77 is a parameter along the particle's 
trajectory as defined by the qap metric, z^ are the 4 
functions describing the particle's position on the brane, 
and Mp is the particle's mass parameter. Using (2.13 we 
find the stress-energy tensor 



(54(Z^ - z^) 



dzP dz^ 

=-^qapqf3\-r^-r^d7j 

q di] drj 



(7.2) 



The contribution from the particle to the total action is 
1 



^±_Mp 



, az^az^ 
'^"^ dr, dr,'^''- 



(7.3) 



Varying this with respect to the trajectory Zp and de- 
manding that rj is an affine parameter yields that the 



rr.± 



f 



C±Ep- 



;UaUf3S{p - Pp) 



1 

47r 



00 / 

EE 

l—l rn— — l 



(7.7) 



Here, we have defined 
{GMf 



Mpe'^y± ' 



9a0u"^^^y e(")=5t. (7.8) 



As usual, E is the particle's energy per unit rest mass 
defined with respect to the timelike Killing vector . 



Comparing (6.3) and (|6.7c| with (|7.7|), we see that 
/ 



Ai 

A2 



f 



inEp' 
Epp 



UaU[3S[p- Pp{t)], 

2 S[p- Pp{t)], 



A, = 



47r/p2 



inEp- 



S[p- Pp{t)], 
J[P- Ppi^)], 



(7.9a) 
(7.9b) 
(7.9c) 
(7.9d) 



where pp — dpp/dr. Here, we have identified L as the 
total angular momentum of the particle (per unit rest 
mass), defined by 



I? 



a B 
la/3U U*^, 



L 



L 
GM' 



(7.10) 
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Note that for particles travelling on geodesies, E and L 
are constants of the motion. These are commonly re- 
parameterized [H] in terms of the eccentricity e and the 
semi-latus rectum p, both of which are non-negative di- 
mensionless numbers: 



(p-2-2e)(p-2 + 2e) 



P{P 

2 



3-e2) 



(7.11) 



The orbit can then be conveniently described by the al- 
ternative radial coordinate x, which is defined by 



P = 



P 



1 + e cos X 



(7.12) 



Taking the plane of motion to be = tt/2, we obtain two 
first order differential equations governing the trajectory 



dx 
dr 

d(j) 
d^ 



(p — 2 — 2ecos x)^(p — 6 — 2ecosx) 



-,1/2 



p4(p_2-2e)(p-2 + 2e) 



p{p — 2 ~ 2ecos 
p4(p_2-2e)(p-2 + 2e) 



-,1/2 



(7.13) 



These are well-behaved thorough turning points of the 
trajectory dpp/dt — 0. When e < 1 we have bound orbits 
such that p/{l + e) < pp < p/{l — e), while for e > 1 we 
have unbound 'fiy-by' orbits whose closest approach is 
Pp = P/{^ + e). To obtain orbits that cross the future 
event horizon of the black string, one needs to apply a 
Wick rotation to the eccentricity e ie and make the 
replacement x *X + ■'i'/2. Then a radially infalling 
particle corresponds to e = oo. 

Since this type of brane matter will be the topic of 
the rest of this paper, it is worthwhile to write out the 
associated source terms in the wave equation explicitly 
as a function of orbital parameters 



'2pPp 



3V47r£;(2 + p?p^) 
6pE'^ f /iV^ - 4 



/ VmV + 2 

PPp „ 



{2p^ + ?,L^)5'[p-pp{T) 
6[p-Pp{T)] , 



5. 



Note that 



Q^J^-nEp 

\pp\ < L 
Pp = ^ 



S'[p-Ppir)]. 



5^ = 0, 



(7.14) 



(7.15) 



S > 1 ^ 5^ > 5, 



That is, the particle's speed is always less than unity, the 
sources wave equation vanish if the particle is stationary 
or in a circular orbit, and high-energy trajectories imply 
that the system's dynamics are not too sensitive to brane- 
bending modes ip ip. 



VIII. SOME TYPICAL WAVEFORMS 

We shall now integrate our coupled system of equa- 
tions (6.14) (i.e. the master equation and the brane bend- 
ing equation) for a variety of different orbital parame- 
ters. Before we investigate some of the typical wave- 
forms which appear, let us briefly digress on some issues 
involved in the integrations and the resulting waveforms. 

First let us consider the solution to the wave equa- 
tion, (6.141, with no source, 5 = 0. This solution is com- 



monly excited by taking Gaussian initial data somewhere 
near the photon sphere and letting the system evolve. 
For a normal black hole this results quasi-normal ring- 
ing followed by a power law tail at late times, as seen by 
a distant observer. In our case, quasi-normal ringing is 
subdominant and instead the signal behaves roughly as 
an oscillating power-law 



sin [u;(r, p)t], 



(8.1) 



where a and uj are slowly- varying functions of r, (as 
compared to the characteristic timescale l/p). While it 
would take a detailed numerical investigation to deter- 
mine the precise nature of these functions for the S-wave 
potential, it has been shown |M] that at late times 
T ^ l//i^ we have, independently of /i, a — > 5/6 (from 
above) while the frequency increases to its asymptotic 
value p as follows: 



a;(r,/i) p 



2/3' 



(8.2) 



The next issue concerns our approximation for the i5- 
functions appearing in the source terms (7.141. These 
we shall approximate by a Gaussian profile in the x- 
coordinate E5l 1261: 



5{r - rp) 



1 



GMf 



6{x — Xp) 
1 

GMfy/^S 



exp 



(x -^p) 



(8.3) 



which becomes exact in the limit e 0. Provided that 
the full width at half maximum (FWHM) 2Vln2£ of the 
Gaussian is much less than the characteristic scales we 
are interested in - namely p and the 'width' of the po- 
tential - this is a good approximation (and is also why 
we choose a Gaussian in x and not r, so that the profile 
remains thin inside the photon sphere). 

The third involves our treatment of initial data. In 
normal relativity, one switches on the interaction at some 
time; the shock in the wave equations produced by this 
propagates way at the speed of light. Here, however, 
a massive mode signal is produced which decays very 
slowly. This makes it difficult to disentangle the real 
signal we are interested in from this spurious signal; we 
shall discuss this further as it arises. 

A further issue which appears is the gravitational 
waves produced just by the unaccelerated motion of the 
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particle itself. Evolving a geodesic compact source in flat 
space within GR does not produce gravitational waves 
(at linear order). With massive modes of the graviton 
present, however, this is not the case: an observer sees a 
wavetail after the particle has passed with a wavelength 
roughly that of the massive mode. This effect also tan- 
gles itself up in the waveforms we are actually interested 
in. 

These issues are illustrated in Fig. |2] where we show V' 
for a particle moving in the far field as seen by an observer 
located at a; = 200. The particle is on a plunge orbit with 
e = 0.1 and p — 3.09, and we have shown the lowest mass 
mode, mi — 0.5. The integration was started with the 
particle located at a; « 240, and we have used e — 1/3. 
There are three key features. The first is the bow wave 



particle 

location 



amplitude 



time 



bow wave 
+ 
junk 



massive mode tail 




» particle's 
\J\ Gaussian 




FIG. 2: Waveform of a particle moving past an observer in 
the far field. 

which precedes the particle: this is just junk from the 
initial data which we want to minimise. This junk will 
not interfere with the signal from the particle interacting 
with the black string, provided we start the simulation 
when the particle is in the far field: in this case, the spike 
from the particle increases by over an order of magnitude 
by the time it gets to x ~ 30, dwarfing any contamination 
from the junk. 

The second feature is the particle passing the observer: 
the disturbance may be compared to the width of the 
Gaussian, whose FWHM is displayed by the width of the 
stem of the arrow pointing to it (and is thinner than 
the width of the line displaying the signal). We can see 
that the disturbance length scale is much wider implying 
that the Gaussian is thin enough. The main part of the 
signal is the massive mode tail which exists in the par- 
ticle's wake. This has a characteristic power law decay 
discussed above; this part of the signal causes problems 
later as typically it will not have decayed away by the 
time the signal from the black string reaches the observer, 
for interesting observer locations (we shall see that inter- 
esting signals occur relatively near the black string, so 



this has ramifications later). 



Plunge orbits and the hierarchy of massive 
modes 



We shall investigate here in some detail the situation 
of a plunge orbit depicted in Fig. [3] with e — .\ and 
p = 3.09 (which corresponds to £: = 2.0, L = 9.78). We 



particle path 



crosses photon 
sphere at t= 70 



particle at x=50 



FIG. 3: Schematic of our plunge orbit. 

shall show how the hierarchy of mass modes contribute to 
the total spherical signal. Assuming we are in the "light 
mode" regime m„£ ^ 1, the KK masses are given by 



Mr, 



[1,1.831,2.655,3.477, 



(8.4) 



The string of KK masses we shall use has /ii = 0.5, which 
starts just above the GL instability, where ^.ql ~ 0.4301. 
For d/l = 20 this corresponds to a black string of mass 
4.3Mo, while for d/i = 35, we have a 1.4 x IO^Mq black 
string. We will present composite solutions for ip; i.e.. 



(z+) i)n{t,x), 



.5) 



where tpn is the numeric solution for a given mass //„ and 
the z+ parameters are given by (4.151. We have assumed 



that both the observer and the source are on the visible 
brane. Note that if we wanted to reconstruct the full 
spherical GW sign al, we would first have to apply the 

(n,s) 



inversion formulae (6.191 to each of the to get h)^ 



[c.f. (6.8)] and then sum over n using (6.2 1 to obtain the 
spherical part of hAB- However, the simplified compos- 
ite signal given above will capture most of the essential 
features of the complete spherical GW signal, and will be 
sufficient for the qualitative discussion given here. 

We show, in Fig. [4] the composite signal ill, and the 
brane-bending contribution ^, as seen by an observer at 
X = 100 for this plunge orbit, starting when the particle 
passes the observer at r = 0. The integration was started 
with the particle located at x « 240, so initial data prob- 
lems give a very small contamination of the signal, and 
we have used e = 1/3. 

The gravity wave signal, -0, has two distinct parts. The 
first, as we have discussed, is from the particle itself and 
the wave-tail it leaves in its wake. But now we have con- 
tributions from higher mass modes, which give a distinct 
wobble to the tail, shown in the bottom left blow-up. The 
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black string signal 




^ , ^ , ^ , > 

200 400 600 800 1000 



FIG. 4: Composite signal from a plunge orbit (top) with brane 
bending contribution (bottom) for an observer at a; = 100, for 
the first 4 mass modes. 



second is from the black string itself (upper blow-up) . As 
the particle falls into the black string it emits radiation 
which in the frequency domain is sharply peaked about 
the frequency of the massive mode (of which more later) . 
Being massive, much of this radiation falls into the string, 
but some of it makes it out to the observer, the first hint 
of which arrives around r « 300. This signal reaches a 
peak around r ~ 4 — 500 - a considerable length of time 
compared to a comparable GR signal - and then gen- 
tly turns into the characteristic power-law tail fall-off. A 
key feature of this is the lack of influence the higher mass 
modes have, compared to /ii; the signal with /ii removed 
is shown in the blow-up at bottom right. We can see 
from the relative scales that this is suppressed by nearly 
two orders of magnitude. Compare this to the earlier 
tail from when the particle passes the observer - mas- 
sive modes higher than /ii are clearly visible there. The 
conclusion of this is that signals arising from excitations 
of the string itself are overwhelmingly dominated by the 
lowest mass mode. 

The overall amplitude of the excitation is worth not- 
ing: ^inax '--^ 5 X 10^*. Given that the source term from 
the particle is 0(1), one might naively expect a signal of 
comparable strength - indeed, this is roughly what hap- 
pens in GR. Such a weak excitation clearly indicates that 
spherical massive modes are only weakly stimulated by 
the particle in-fall. In part this is due to the fact that 
some of the signal falls into the string; more on this later. 

Finally, we come to the composite brane-bending con- 
tribution to the signal. The signal, which is indepen- 
dent of ^, is pretty featureless. As the particle passes, 
a dent in the brane accompanies it; this reaches a peak 
after the particle has passed, and slowly relaxes back to 
zero without oscillating. As the brane remains signifi- 
cantly bent long after the particle passes, this extends 



the total source feeding the gravity wave signal beyond 
the particle's Gaussian. Thus, the black string gets a 
far longer stimulation than it would otherwise get from a 
point source: the brane bending signal is partly responsi- 
ble for the length of time ip remains peaked in the latter 
part of the signal. This may be seen by the fact that the 
tail part of the signal has a power-law fall off of a ~ 1.1 at 
T K, 1000, so hasn't yet reached the asymptotic late time 
value of 5/6. However, comparison of the signal with the 
brane bending switched off reveals that the contribution 
to the amplitude is only of the order of a few percent. 

B. Bound flower-shaped orbits: steady state 
waveforms 

Let us now investigate the signal which comes from a 
bound orbit, illustrated in Fig. [5] and explore how the 
signal changes with distance from the black string. We 
concluded from the previous section that the higher mass 
modes add only a small correction to the full signal, so 
here we shall only investigate the signal arising from the 
lowest mass, /ii — 0.5. We choose an orbit with e = 
0.5, p — 7.05, and we take e = 1. 

periastron 
x=5.3 




FIG. 5: Illustration of the flower-shaped bound orbit used for 
investigating the steady state waveforms 

Integration of the equations is complicated by initial 
data, once again, but more-so than in the plunge case. 
This is because the amplitude of the source in the wave 
equations increases with decreasing x, so the best we 
can do is start the integration at the apastron where it's 
smallest, and wait for the contamination to pass the ob- 
server. Unfortunately the wave tail makes this quite a 
long time - roughly r ^ 4 — 5000 for an observer around 
X < 100, compared to r ~ 100 in GR. After this time 
the desired steady-state waveform is reached, which we 
set to r = 0. 

In Fig. [6] we show the results of this integration for 
several observer locations. As the source evolves it oscil- 
lates along the x direction, between about 5.3 and 17.7. 
It can be immediately seen that observers see radically 
differing signals depending on whether they are in the 
near, intermediate, or far zones. 

Near Zone: An observer sitting near the photon sphere 
around x = Q will see a relatively normal signal: 
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a gravity wave propagating at light speed (since 
y « there) which precisely mimics the behaviour 
of the source. Around x = 10 - located in the 
'middle' of the orbit - the amplitude of reflects 
the source passing back and forth. 

Intermediate Zone: Here things get a bit more inter- 
esting. As we move further out to a; '--^ 40 the direct 
orbital signal gets damped dramatically - by 4 or- 
ders of magnitude. This is due to low frequency 
parts of the signal being exponentially damped as 
they travel through the potential (all frequencies 
roughly less than ji get damped). By a; ~ 30 the 
massive modes become visible as a high frequency 
wobble riding on the orbital part. And around 
x ~ 40 the two components become equally domi- 
nant. 

Far Zone: As we observe from more distant locations 
where the potential is almost flat virtually all the 
low frequency components have been suppressed, 
and we are left with a low amplitude massive sig- 
nal. A gentle oscillation to the envelope is all that 
remans of the orbital signal. 

One of the interesting and unexpected things about 
this is that the massive modes are excited at all given 
that the frequency of the source is orders of magnitudes 
smaller than the mass. Also note that the amplitude of 
the signal falling into the black string is about five orders 
of magnitude larger than the signal which makes it out. 



C. Flyby orbits 

Finally we shall consider the case of unbound orbits 
with high angular momentum. We choose an orbit with 
e = 3.0 and p = 12.00001 (see Fig. fig:fly). This com- 
pletes about four orbits of the black string, very nearly 
touching the photon sphere at is periastron (t = in this 
simulation). If we think of this orbit in terms of the x 




FIG. 7: Schematic of the fly-by orbit. 

coordinate, the particle comes in on a nearly null ray, de- 
celerates very rapidly to zero at the photon sphere where 
it lingers for r ~ 150. It then rapidly accelerates off to in- 
finity (of course the particle is geodesic so experiences no 
acceleration). The source in the wave equation becomes 
very large at these two events. We see in Fig. [8] how this 
acceleration can induce a strong GW signal. First con- 
sider the GW which falls into the string. This consists of 
two pulses corresponding to these two accelerations. As 
the potential is basically flat in this zone the pulses re- 
flect closely the time evolution of the source. We can see 
some massive modes on top of this caused by a reflection 
off the potential. 

In the far zone, on the other hand, only the second 
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particle passes observer 



>- black 
string 
signal 




time 



FIG. 8: The GW signal from a fly- by orbit, for the lowest 
mass mode. 



acceleration produces a GW signal - and this is of nearly 
the same amplitude as the signal which falls into the 
string. In the previous two examples the far zone signal 
was orders of magnitude smaller than that which passes 
the horizon. The burst of GW seen at a; = 200 nearly 
coincides with the passing of the particle, although the 
blow-up shows the the particle precedes the signal some- 
what. The peak part of the signal has some interesting 
wobbles; after this the waveform takes its familiar shape 
of a decaying massive mode. 

The interesting aspect of this simulation is that the 
amplitude is orders of magnitude larger than the previ- 
ous two cases. We can understand why this happens by 
examining the source term, given by Eqs. (7.141. The 



two terms which are most important are the coefficients 
of 5'[p — Pp(t)]. This is proportional to vT , which when 
E and L are large scale like {p — i ~ e^) 
the simulation above); thus when p — > 3 - 
infinite source and signal. Note that u'' = at the peri 



astron, so the source peaks near the periastron, when r 
is small, but not far away when 1/r terms kick in. We 
have numerically confirmed that the signal which makes 
it to infinity does in fact scale in this way. 



IX. DISCUSSION 

In this paper, we have presented the derivation and 
numeric solution of the equations of motion for gravi- 
tational waves in the black string braneworld sourced by 
brane localized matter. In j jlTj we presented a generalized 
two brane Randall-Sundrum model and then specialized 
to the black string background. In §III| we considered 
the linear perturbations of the model and the introduced 
the Kaluza-Klein massive mode decomposition pV| The 
limit under which the model reduces to Brans-Dicke the- 
ory was discussed in SjV] which led to a constraint on 
the brane separation. We discussed the specialization of 



1/2 (= 102 for 
e2 we have an 



the formalism to spherical radiation (^VIl and pointlike 
sources (j |VII I. Finally, in ^ VIII we presented the results 
of numeric simulations of the spherical GWs produced by 
perturbing bodies undergoing plunge, bound, and fly-by 
orbits. 

Future work on this model involves improving our sim- 
ulation techniques by incorporating characteristic inte- 
gration techniques [8] that more naturally deal with the 
delta-functions in the GW source. It would also be in- 
teresting to consider more realistic modeling of sources 
of finite size. Once this is accomplished we can build up 
a bank of simulations for a variety of orbital parameters, 
choices of /i, and other multipoles. One can then sys- 
tematically begin looking for these waveform in the data 
obtained from gravitational wave detectors, and thereby 
provide a means of further constraining the Randall- 
Sundrum braneworld model. 
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